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Optimal Index Assignment for 
Multiple Description Scalar Quantization 

Guoqiang Zhang, Janusz Klejsa, and W. Bastiaan Kleijn 

Abstract 

We provide a method for designing an optimal index assignment for scalar if-description coding. The method 
stems from a construction of translated scalar lattices, which provides a performance advantage by exploiting a 
so-called staggered gain. Interestingly, generation of the optimal index assignment is based on a lattice in K — 1 
dimensional space. The use of the K — 1 dimensional lattice facilitates analytic insight into the performance and 
eliminates the need for a greedy optimization of the index assignment. It is shown that that the optimal index 
assignment is not unique. This is illustrated for the two-description case, where a periodic index assignment is 
selected from possible optimal assignments and described in detail. The new index assignment is applied to design 
of a if-description quantizer, which is found to outperform a reference ^-description quantizer at high rates. The 
performance advantage due to the staggered gain increases with increasing redundancy among the descriptions. 

Index Terms 

Multiple description quantization, Index assignment. 

I. Introduction 

THE real-time transmission of multimedia content over contemporary packet-switched networks generally 
requires a coding scheme that can address the effects of packet loss. Multiple-description coding (MDC) 
creates a plurality of descriptions of a source signal. The descriptions are embedded into packets and transmitted 
over a lossy network. The descriptions are mutually refinable and any subset of them can be used to reconstruct the 
source signal. The more descriptions reach the destination the better the reconstruction is. This property naturally 
requires the introduction of redundancy, which is used to adjust the trade-off among the distortions arising from 
the possible description loss scenarios. 
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Multiple description quantizers are the most prominent class of MDC schemes, since they are almost directly 
applicable in many practically relevant source coding scenarios. The first multiple description quantization schemes 
were proposed by Vaishampayan for a two-description scalar case CO, El- The scalar schemes were extended to the 
two-description vector case in 0, El, 0. A symmetric A-description vector scheme was proposed by 0stergaard 
et al. in 01 and extended to even more general asymmetric case in Q. 

A majority of MDC quantization schemes is based on the construction of a fine "central" quantizer and A coarser 
side quantizers. The central quantizer is related to the side quantizers by means of an index assignment mapping, 
which uniquely maps a quantization index to a tuple of A indices that are associated with the side quantizers. The 
design of such an index assignment plays a central role in optimizing the performance of MDC. 

The research on index assignments for MDC is vast. The first practical index assignment schemes were proposed 
in HI- The design of an efficient index assignment can be formulated as a combinatorial problem (H, [9), a graph 
optimization ifTOl , a transporation problem in operations research ifTTTl . If lattice codebooks are used, the design 
problem of the index assignment can be formulated as a labeling problem, where each lattice point of the central 
quantizer is mapped to a unique if -tuple consisting of points belonging to A respective lattice side quantizers. The 
labeling problem was solved in (6l, Q by exploiting properties of clean sublattices and using bipartite matching. 
An efficient index assignment scheme was proposed in lfT2ll . where an auxiliary so-called fractional lattice was 
introduced significantly simplifying the labeling problem. The fractional lattice was also used in [11], where an 
index assignment based on a transportation model was developed. 

This paper aims at designing an optimal index assignment for if -channel multiple-description scalar quantization. 
The proposed index assignment scheme utilizes a so-called staggered gain, which stems from a construction of 
translated lattices. Several practical multiple description schemes attempt to utilize this gain in a heuristic manner 
for a two-description scalar case lfT3l . a two-description vector case lfl4l and A'-description case with a two-stage 
coding lfl5l . In our work the index assignment exploits the staggered gain for an arbitrary number of descriptions. 
The staggered gain becomes considerable in the case of high redundancy among the descriptions. However, the gain 
vanishes as the redundancy decreases and the performance of the proposed index assignment becomes equivalent 
to existing schemes (e.g. lfT2lO . Nevertheless, the new index assignment is generally advantageous, since the cases, 
where the redundancy vanishes, are of low practical importance. 

In the proposed scheme, a central quantizer and A side quantizers are designed to be translated Z 1 lattices. A 
simple mapping function from a central point to a A-tuple of the cartesian product of the A side quantizers is 
proposed. The selection of a good A'-tuple to label a central point is performed by choosing a point of a translated 
Ak-i lattice with a short distance from the origin. The index assignment is shown to be optimal under a common 
decoding process. 
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In this work, we use a so-called reference quantizer that resembles the idea of a fractional lattice |[T2l . which 
facilitates the design of the index assignment. A difference to work of 1 12 ] and iTTTTl is that the reference quantizer 
arises from a construction of translated side quantizers, which utilizes the staggered gain. We show, that using 
the proposed reference quantizer, we only need to label the central points in a reference quantizer cell. Further, 
the labeling operation becomes straightforward by using the Ak-\ lattice. Thus, the labeling complexity for the 
proposed index assignment is reduced compared to that of |[T2l and iTTTTl . In addition, the performance for the case 
of non-vanishing redundancy is improved. 

The use of the Ar-i lattice in designing of the index assignment has a number of practical advantages. It 
facilitates generation of the index assignment at hand, leads to low operational complexity of an MDC scheme and 
provides analytic insight into the performance of the obtained multiple description quantizer. We illustrate these 
properties by describing a regular index assignment for the two-description case in more detail. The index assignment 
based on the Ak-\ lattice is applied to the i^T-description scalar quantization. The performance advantage over 
|[T2l is demonstrated. The high-rate performance is further evaluated by considering a distortion product |[T6l . 

II. Preliminaries 

Suppose a source random variable V is to be encoded and transmitted through K channels. Denote its realization 
as v. It is first quantized to the nearest point A c = Q(v) of a central quantizer A c - K descriptions of A c are then 
produced and transmitted through K separate channels. Assuming symmetry of the channel conditions, we consider 
designing balanced descriptions, where the transmission rate is the same per channel and the decoding distortion 
only depends on the number of received descriptions. The mean squared error (MSE) is taken as the distortion 
measure. 

Each description is a quantization index describing the associated side quantizer point. We denote the K side 
quantizers as Ai, i = 0, . . . , K — 1. The descriptions are produced through an injective labeling function a : A c — > 
Ao x A\ . . . x Ak-i, expressed as 

a(A c ) = (A , Ai, . . . , Xk-i), h € A- (1) 

We write each component function of a as ati, i = 0, . . . , K — 1. The labeling function usually results in that the 
central cells associated with each side quantizer point is disjoint. Upon receiving all the descriptions, the central 
point A c is determined uniquely by the inverse mapping a -1 . This requires that each i^-tuple (Ao, Ai, . . . , Xr-i) 
is used at most once. In principle, there are 2 K — 1 decoders as there are that many possible channel states. As 
the channels are symmetric, we consider designing a MDC scalar quantizer such that the decoding operation and 
the (mean) distortions are only affected by the number of received descriptions. To achieve this goal, we exploit 
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a common decoding process Suppose k (1 < k < K) out of K descriptions are received. Considering which 
k description are received, there are (^) different configurations. Let C^ K,K ^ denote the set consisting of all the 
possible configurations. We denote the K-tuple associated with an element I € C^ K,K ' as {Aj., j = 1, . . . , k}. The 
reconstruction of the source V for some I € £(^> K ) is taken as the average of the received descriptions |6): 

V = ~it X lr ( 2 ) 

i=i 

Strictly speaking, the estimator V in (O might not be optimal. To achieve optimal estimation, the central cells that 
contribute to the element I must be known, complicating the design problem. Conversely, the use of (O provides 
a good estimate and facilitates the design of index assignment. Under the averaging operation in (0, the decoding 
system is essentially simplified to two decoders, see Fig. [T] The quantity denotes the distortion when k 

descriptions are received. D^ K K ^ is referred to as a central distortion, and > K = 1, ■ ■ ■ , K — 1 as side 

distortions. 

\K,K) 

Fig. 1. The schematic MDC scheme of the proposed index assignment. The quantity n indicates the number of received 
descriptions. 

Upon specifying the decoding process, the challenge is to design an index assignment (specify the labeling 
function a) to achieve good performance. Optimality and simplicity of the index assignment are the two main 
concerns in our work. 

III. Index Assignment 

We first describe the arrangement of a central and K side quantizers. Then we formulate optimality of an index 
assignment. After that, we present the proposed optimal index assignment. We then study the properties of the 
index assignment. Finally, we consider the index assignment for the two-channel case in the form of an IA matrix. 

A. Setup of Central and Side Quantizers 

The side quantizers are constrained to be translated Z 1 lattices. For the case of K > 1 descriptions, we define 
the side quantizers as 

Ai = {K( ■ Xi + (2i - K + l)^\ Xi £Z}, i = 0,...,K - 1, (3) 

where the scaling parameter £ is introduced to adjust the distances between the side quantizer points. We refer 
to Xi in © as the coordinate of A%. The K side quantizers are arranged so that Ai+\ is obtained by translating 




k < K 



iveraging 



Ai leftward by of the (side) cell width. Intuitively speaking, the arrangement produces a joint quantizer that 
achieves lower distortion than that of each side quantizer. The performance improvement due to translating side 
quantizer to produce a finer joint quantizer is referred to as a staggering gain lfl~7l . See also [18] and I19| for MDC 
schemes exploiting translated lattices. 

The joint quantizer obtained by combining the K side quantizers is referred to as the reference quantizer, denoted 
as A r . The reason for introducing the reference quantizer is to separate the arrangement of a central quantizer from 
the K side quantizers. The centroids of the reference cells take the form 

A r = {(- z\z G Z}. (4) 

It can be easily shown that 



Ar — ^ ^ Ai, (5) 

i=0 

which states that the centroid of any .fT-tuple is a point of Ar- We point out that the fractional lattice introduced in 
lfT2l has the same property as (f5]). Specifically, in lTT2l . the fractional lattice is defined to be a side lattice scaled by 
1/K, and without lattice translation. In our work, the reference quantizer arises naturally from the arrangement of 
side lattice quantizers. The reference cell width is £ as compared to the cell width K( of a side quantizer. Thus, 
each side quantizer cell contains K points of the reference quantizer. 

T H i ' t i ' 8 i ' 1 i ' 1 

reference quantizer 
d 




central quantizer 



Fig. 2. The quantization structure for K = 3. The points denoted by A, V and represent the points of the three side 
quantizers. 

Based on the reference quantizer, the central quantizer is defined as 

A c = {^-y + ^-mod(M + l,2):yeZ}, (6) 

where M is an integer. Thus, the cell width of the central quantizer is -j^. The definition guarantees that there are 
M central points within each reference quantizer cell and also that the distribution of the points within each cell is 
the same. From © and (J6J), it is immediate that there are KM central points within each side quantizer cell. We 
define N = KM as the redundancy index. The parameter M (or N) serves as a trade-off factor between central 
distortion and side distortion, which can take any value of Z + . An example of the quantization structure for K = 3 
is illustrated in Fig. [H Note that the K side quantizers and the central quantizer are arranged periodically along 
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the line R. We only have to label the central points in a reference quantizer cell, which we will discuss in next 
subsection. 

The Voronoi region V(X) of a point A in a quantizer A is defined to be 

V(X) = {x\(x — A) 2 < (x — A) 2 , VA G A}, (7) 

where the ties are broken in a systematic manner. To describe the relation between the reference and central 
quantizers, we define a discrete Voronoi region associated with each A r G A r as 

Vr(X r ) = {X c G A C \(X C - X r ) 2 < (A c - A^) 2 , VA' r G A r }. (8) 

The definitions of A c and A r guarantee that the central points do not lie on the boundaries of the reference cells. 
From ©, the cardinality of V r (X r ) is |T^(A r )| = M for any A r G A r . 

By studying Q, the points of the side quantizers are periodically and evenly distributed over the line M. Now we 
formulate this periodicity w.r.t. the reference quantizer points. We assign a quantizer-tuple (Ao,Ai, ■ ■ . , Ak-i) to 
the element G Ar to describe the geometrical relationship between the K side quantizer and A r = 0. Considering 
the periodicity of side quantizers, the quantizer-tuple for a reference quantizer point A r = £ • z is 

7(C ' z ) = {•^■{mod(z,K)}, ■ ■ • ) •^■{mod(z+K~l,K)}) ■ (9) 

The 7(-) function in © exhibits periodicity along with reference points £ • z with period K. This geometrical 
property facilitates the design of optimal index assignment and further the generation of the balanced descriptions. 
We now proceed with a definition of optimality of an index assignment. We first define a cost 



B 



(X c ,k)= (A c --^A /j ) 2 (10) 



K 

; e £(jr,«) j=o 



to evaluate a labeling function for a particular central point A c for the case of k descriptions received. Note that 
the defined cost is a geometrical measurement, and is unrelated to channel conditions. The following theorem 
decomposes the cost (fTOl (see lTT2l . |20l for details). The result applies to any dimensionality as long as the l 2 
norm is taken to measure the error. 

Theorem 3.1: |[T2l Suppose A c G M 1, is associated with a iC-tuple (Ao, Ai, . . . , Xk-i), where Aj G M L , i = 
0,...,K-1. Then for any 1 < k < K: 



-riK.K) o — 1 v ' 



A, - A II 2 



K-l 

(11) 



v ; »=0 



where || • || 2 denotes the I2 norm and A denotes the centroid of the -fT-tuple, i.e. A = T^f=Q ^j- 

The two terms on the right side of (fTTT) can be interpreted geometrically. The first term measures the squared 
distance (SD) between A c and the centroid of a K-tuple. The second term computes the sum of squared distances 
(SSD) between the components and the centroid of the if-tuple, which captures the geometrical structure of the 
K-tuple itself. 

The optimal index assignment for the case of k (k < K) descriptions received is defined as 

a^=rmn-±- £ B(X c ,k), (12) 

where \A C \ represents the cardinality of the central quantizer. We note that B(X C , n) is a function of k. By applying 
(fTTT) . (fT2l can be rewritten as 



(K,k) 



A'-l 

(13) 



Kk(K 

v ' A c eA »=o 

where A(A C ) = X^o* a «(Ac)- The right-hand side of equation ([T"3l consists of two parts, the SD and SSD 
costs. We propose an index assignment that is optimal with regard to any k < K in (fT"3l) . We show that the index 
assignment minimizes both the SD and SSD costs simultaneously. 

B. Construction of Optimal Labeling Function 

Similarly to the two-channel case [lj, the index assignment for the general ET-channel case can be posed as 
a geometrical minimization problem. We show that once the central points in V r (0) are labeled, the rest can be 
obtained easily by exploring the regularity of the central and side quantizers. The index assignment is then shown 
to be optimal w.r.t. (fT3"T ). 

From ©, it is known that the centroid of any i^-tuple used for index assignment is a reference point. Thus, all 
the if-tuples can be grouped with respect to their centroids or, equivalently, their reference points. Each reference 
point is associated with many i<C-tuples. We evaluate those K-tuples with a common centroid by measuring their 
SSD costs. We start with a If -tuple with centroid A r = 0. From ([3]), any point of a side quantizer Ai is determined 
by its coordinate X{. Thus, the SSD cost takes the form 

K-l K-l 

Kt -XA + (2i- K+V, 



J2 X i = ^(KC-x. + ^-K + l)^) 2 (14) 



1=0 8=0 

K-l 

K 



1 K ~ l 
subject to — 22 



i=0 



0. (15) 



The constraint is imposed to reflect the centroid condition. Denote X = [ X q x\ ... xk-i ] T > representing a 
coordinate vector of a i^-tuple. Eq. (fT4b-([T5T) can be further simplified as 

J{X) = K 2 C 2 || X + s || 2 (16) 
K-l 

subject to Xi = 0, (17) 

i=o 

where 

s = 7^[-(K-l) -K + 3 ... K-l] T - (18) 



Take X as a point in if-dimensional space. The constraint ( fTTT ) forces the points to be in a hyper-plane of 
dimensionality K — 1. As each component of X only takes integer values, (fTTT ) defines an Ajc-i integral (the 
inner product of any two lattice vectors is an integer) lattice |[2ll . Thus, the cost J(X) in (fT6l ) can be interpreted as 
measuring the squared distance of a point of a translated Ak-\ lattice from the origin regardless of the multiplying 
factor K 2 C, 2 . The translated lattice takes the form 

Ak-i + s, (19) 

where the vector s is the translation. We denote the translated lattice as Ax-i. Thus, one can associate all K-tuple 
candidates that have common centroid A r = with the points of Ak-x- The translated lattice Ak-\ reveals the 
geometrical relation between the side quantizers and € A r . Note that each component Xi of X is the coordinate 
of the side quantizer Ai. This connects A r = with the quantizer-tuple (Aq,Ai, ■ ■ ■ , Ak-i) through X, as is 
specified by the 7(-) function. 

Next, let us consider the SSD cost of a iC-tuple with a reference point ( • z as its centroid, and further the 
geometrical relation between £ • z and the side quantizers. Similarly to (fT4l)-(fT5T), the SSD cost can be expressed as 

K-l K-l t 

J2(\ l -(.z) 2 = J2(K( ■ Xi + (2i - K + - ( ■ z) 2 (20) 



k=0 i=0 
K-l 

K 



1 K ~ 1 
subject to — 



K(- Xl + (2i- K + l)^ 



C • z. (21) 



i=0 

The above two expressions can also be simplified and characterized by the translated lattice A^-i, which is 
described in the following Proposition. 



Proposition 3.2: The SSD cost defined in (120I)- (|211) nas a simplified form 

K-l 



J(X) subject to Xi = 0, (22) 



i=0 
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where 



x mod(z,K) \-K^ 
x mod(z+l,K) ~ l^jr\ 



(23) 



™ f I z+K-l I 

The operation [x\ takes the largest integer not exceeding x. 
See Appendix |A] for the proof. 

Note that the ith component of X in (l23l is the coordinate of the side quantizer -A m od(z+i-i,if )■ This again relates 
the reference point £ • z with a quantizer-tuple through X, confirming Q. Eq. (l22l shows that all reference points 
are geometrically equivalent in terms of SSD cost as J(X) always has the same structure. In other words, for every 
K-tuple with centroid A r = 0, there exists a fT-tuple with a centroid A r = £ • z that gives the same SSD cost. 

The search procedure for good index assignment can be performed in two steps. First, the central points within 
the fundamental reference cell, {A c G V r (0)}, are considered and labeled. The K-tuples exploited are enforced to 
have a common centroid A r = 0. From (l23l . the labeling function for any central point A c + m£ S A c , m € Z, 
A c G V r (0), can then be obtained as 

m + i 



a 



{mod(m+i,K)}(^c + Tn() = «i(A c ) + [ J 



(24) 



where i = 0, . . . , K — 1. The operation defined in (1241 ) guarantees that the K-tuples used to label the central 
points within V r (X r = ( ■ z) have a centroid £ • z. This ensures that no if -tuples are reused by performing (l24l . 
It is seen that the extension of the index assignment from V r (0) to V r (X r ) generally involves both translation and 
permutation (changes of subscripts of a labeling function an). The permutation operation exhibits periodicity with 
period K(. Thus, the extension of the index assignment from Ua =0^ ^r-(A r ) to Ua =o X V r (X r +mK(), m € Z, 
only involves translation. 

As described before, in order to label the M central points in V r (0), we first constrain if-tuple candidates to 
have centroid A r = 0. This ensures that every central point contributes the corresponding minimum SD cost to 
(fT3l . The selection of good isT-tuples for labeling can then be done by choosing points of Afir-i- Specifically, we 
order the points of Ak-i according to their distances from origin. Denote the coordinate vector of the ith point 
as The first M coordinate vectors {Xi, i = 0, . . . , M — 1} are then selected, which give the lowest M SSD 

costs under the centroid constraint. We then define a bijective mapping 



P:V r (Q)^{Xi,i = 0,...,M-l} 



(25) 



to relate the i^-tuples to the central points. Upon selecting a (3 function, an index assignment is fully determined. 
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The cost in ( fT3l ) incurred by the index assignment is given as 



1 



K — K 



M-l 



M 



E ^ + 



Kk{K - 1) 



E J ^ 



(26) 



A c GK.(0) 



i=0 



We now investigate the optimality of the index assignment. 



Theorem 3.3: For the quantization structure defined by ([3]) and Q and the averaging decoding strategy specified 
by ©, the index assignment specified by a (3 function in d25l ) and (1241) is optimal w.r.t. ( fT3l ) for any k < K. 

Proof: The proof is trivial. First, the proposed index assignment guarantees that the SD part of the costs in (fT3l ) 
is minimized. Second, (1241) and d25l ) imply that all the exploited if-tuples have SSD costs not exceeding J(Xm-i)- 
If unused K-tuples are taken to replace exploited ones, the part of SSD costs in ( fT"3l ) would obviously increase. On 
the other hand, the part of SD costs would either increase or remain the same by switching the exploited K-tuples 
of different central points. This shows that the proposed index assignment is optimal for any k < K. ■ 

Note that a different (3 function in (l25l) does not affect d26l ). This implies that there exist more than one optimal 
index assignment when M > 1. The simplicity of the proposed labeling function is due to the proper arrangement 
of the central and side quantizers. 

C. Properties of the Index Assignment 

In this subsection, we first study the distribution of the points of Ak-i- Specifically, its theta series ETI is 
investigated for distortion evaluation. The issue of generating balanced descriptions is then discussed. Finally, we 
consider the possibility of smooth adjustment of the redundancy among the descriptions. 

To analyze the performance of the proposed index assignment, an essential step is to study the properties of SSD 
costs. This motivates us to look into the theta series of Ak-\, which is defined as 11211 



where q = e 7 ™ 2 . It is seen that the theta series captures the SSD costs of all the K-tuples that have a common 
centroid, rendering an infinite series. The theta series describes the distribution of the number of the translated 
lattice points with a common distance from origin. In principle, one can easily access the squared distances of the 
first M points given the theta series, thus determining the part of SSD costs in d26l ). Further, if the theta series 
exhibits regularity (e.g. the number of points can be parameterized by their common distance from origin), the SSD 
costs of the first M points can be computed easily. Let us denote by A* K _ 1 the dual lattice (see Appendix iBl or 
|[2ll for definition) of Ak-i- The translation vector s in (fT8l) is a deep hole j[2TTl of A* K _ 1 . Informally speaking, 
deep holes of a lattice are the points in the space that are farthest away from the lattice points. The theta series of 
Ak-i can be computed easily with the aid of an A* K _ 1 lattice. 




(27) 
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Proposition 3.4: Let A* K _ l (hole) denote the translated lattice that is obtained by translating A* K-1 with one of 
its deep holes. The theta series of Ajc-i is related to that of A* K1 {ho\e) through 

e^_>) = ^e^_ l(h0le) (4 (28) 



See the proof in Appendix |B] 

The theta series of lattices A*, i = 1, 2, 3, at their deep holes are well studied Ell , thus facilitating the analysis 
of Ai, i = 1,2,3. For example, takes the form 

oo 

e Al = 5>^> 2 . (29) 

m=0 

The expressions of the theta series ©a* (hole) (^) for i > 4 remain to be discovered. Note that the theta series reveals 
the information about the number of points with a particular Euclidean distance from the origin. One application 
of the derived theta series is to verify the exploited K-tuples in labeling the central points. The exploitation of the 
theta series in computing the distortions D( K > K \ i = 1, . . . ,K, will be discussed later on. 

Next we examine if the index assignment can produce balanced descriptions. We study the cost incurred by a 
particular set of received descriptions (corresponding to a I G C^ K,K >). Due to the periodicity of the index assignment 
in d24l . we measure the cost 

E ^~lt, X h) 2 (30) 

for a particular / € C^ K,K \ k < K. For the two-channel case, the two components «o and a\ for (J^ = o^-(A r ) 
are related by 

f ai(A c + C) = a (A c ) + C 

I , A c e v r {0). 

[ a (A c + C) = ai(A c ) + C 

This implies that the roles of the two side quantizers in index assignment are changed symmetrically between V r (0) 
and V r (Q. Thus, the two costs corresponding to the two elements in Z^ 2 ' 1 ) are the same, ensuring the generation 
of balanced descriptions. Similar behavior is observed for the three-channel case, which suggests that the index 
assignment for K = 3 also generates balanced descriptions. When K > 3, other techniques are required to produce 
balanced descriptions, e.g. time-sharing coding. 

The parameter M takes integer values. This implies that the adjustment of the redundant information among 
descriptions cannot be performed smoothly. To overcome this issue, one can add an extra central cell with width 
smaller than (,/M within each reference quantization cell lfT3l . The cell width can vary from to C/M depending 
on the redundancy needed. We will not analyze this method here. 
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D. Index Assignment Matrix for the Two-description Case 

For the two-description case, the proposed index assignment can be visualized by transforming it to an IA matrix. 
We provide a method to parameterize the IA matrix. We show that the obtained IA matrix is also optimal under a 
different criterion ll22l . 

To construct an IA matrix from the index assignment, we index the central quantizer points using its coordinate 
y. The coordinates of the central points within V r (0) are 



, M , , M , . M . 

{y = -LyJ, -L T J + 1, • • • , -L T J + m - 1}. 



(31) 



Similarly, we index the points of the two side quantizers (Ai, i = 0, 1) by their coordinates xo and x\. The 
labeling function can then be simplified as a(y) = (xq,xi). The inverse mapping is denoted as a^ 1 (xQ,xi) = y. 
A quantization unit of side quantizer is defined as 

Co{xo) = |J y- 

Thus, a quantization unit Cq{xq) of side quantizer is the union of central quantizer points that map to xq. A 
quantization unit for side quantizer 1 is defined similarly. 
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Fig. 3. (a) The IA matrix for M = 2, where the bandwidth is 4. The elements in the matrix are the central cell coordinates, (b) The 
extracted patterns w.r.t. Side quantizer 0. (c) The extracted patterns w.r.t. side quantizer 1. The number in (■) describes the translation between 
a column(row) pattern and its original quantization unit. 



To determine an IA matrix, we use the following (3 function: 



/% = -LyJ+*') 



Xi, i 



0, 



, M - 1. 



(32) 



In this situation, the redundancy index N = 2M becomes the bandwidth of the IA matrix. An example for M = 2 
is given in Fig. [3] (a). It is seen from the figure that any quantization unit from both the quantizers is in fact 
a translated version of a fundamental pattern {0,1,3,6}. Fig. [3] (b) and (c) display the translated quantization 
units of the two quantizers, separately. Similar properties are observed for other M. As will be shown below, the 
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fundamental patterns for different M can be captured by an expression in terms of M: 

P(M) = {0} U {i(M - 1 U {(M - l) 2 + M} (33) 

U{(M - l) 2 + M + i(M + l)}^ 1 . 

where M > 1. An IA matrix can be built on P(M). 

We describe the scheme of how to extract the expression P(M) from the specified f3 function. For the two- 
channel case, any 2-tuple can also be related to a line segment connecting two side quantizer points. The length 
of a line segment is the associated SSD cost. This implies that any line segment with length belonging to 

{C,3C,...,(2M-1)C} (34) 

plays a role in the index assignment. Without loss of generality, we consider the quantization unit Co(0). There 
are 2M exploited line segments with one common endpoint Ao = — C/2 (equivalently xq = 0). In particular, any 
element in the set in d34l is associated with two exploited line segments, one having the left endpoint —C/2 and the 
other having the right endpoint —C/2- The corresponding 2-tuples are {(—C/2, C/2 + 2Cxi)|xi = — M, . . . , M — 1}, 
where x\ is the coordinates of the points of side quantizer 1. By computing the centroids of these 2-tuples, the 
coordinates of the resulting reference points are related to x\ by z = x\. When z increases from — M to —1, the 
length of the associated line segment decreases, implying 

a _1 (0,i + 1) - a _1 (0,i) = M — 1, i = -M, . . . , —2. 

The lengths of the line segments when z = — 1 and z = are equal to C- This implies that a _1 (0, 0) — a _1 (0, — 1) = 
M. Conversely, when z increases from to M — 1, the length of the associated line segment increases, implying 

a~ 1 (0,i + 1) - a _1 (0,i) = M + 1, i = 0, . . . ,M - 2. 

Combining the results from these three situations produces the pattern P(M) shown in (l34l . Due to the symmetry 
of the side quantization points, other quantization units Cj(xj), i = 0, 1, give the same pattern. 

We now present how to construct the IA matrix from P(M). Without loss of generality, assuming Co(0) is 
expressed by the pattern P(M) (usually, Co(0) is translated from P{M)). Then Cq{xq) is described by 

P(M) + 2Mx . (35) 

The quantization units are arranged column-wise along the principle diagonal in a matrix, as shown in Fig. [3] (a). 
The constructed matrix implicitly determines the quantization units of side quantizer 1, which can be described as 
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Ci(xi) : P(M) + 2Mxi + M. 

The IA matrix is systematically parameterized with the aid of d35l ). This facilitates the generation of the IA 
matrix. Due to the simplicity of P(M), one can analytically investigate its asymptotic performance (e.g. ll23lD . or 
compute the side distortions and further derive optimal bandwidth in response to varying channel conditions. 

An alternative principle of designing an IA matrix is to minimize the so-called spread, the difference between 
the minimum and the maximum central indices within each quantization unit, subject to a constant bandwidth 12. 
A lower bound on the spread is derived in |[22ll . given as bib — l)/2 where b is the bandwidth. It is immediate from 
d34t that each quantization unit achieves this bound. It can also be shown that the IA matrix produced from any 
other p function also achieves this bound. Thus, the derived IA matrix is also optimal from a viewpoint of spread 
measurement. 

IV. Evaluation of the Index Assignment 

In this section we evaluate the proposed index assignment. For brevity, we name the index assignment as Ak~\- 
based IA. We mainly focus on the high-redundancy case, or equivalently, the index assignment with small M. This 
is because the high-redundancy case is more relevant in practice. 

We consider encoding the Gaussian source X ~ iV(0, 1), as an example. In the past, many theoretical results 
have been obtained for the Gaussian source. Thus, by choosing the Gaussian source, we are able to study the 
performance loss of the ^4/^_i-based IA. 

One popular way to evaluate an index assignment is to analyze its performance under the high-rate assumption. 
By doing this, analytic expressions can often be obtained for approximating the corresponding side and central 
distortions. We point out that in most practical situations, both the irregularity of the index assignment and the 
high-rate approximation negatively affect the accuracy of the side distortion expressions (see 0, lfl2l . 11241 . 0). It 
may happen that different index assignments may have the same expressions for the side distortions, even though 
their real performance is different. 

In the following, we show that the ^i-based IA (i.e., K = 2) is regular. In other words, the accuracy of the side 
distortion expression is only affected by the high-rate approximation. This is due to the fact that the theta series of 
A\ can be nicely parameterized (see d29ll). On the other hand, the theta series of A* K _ V K > 2, take complicated 
forms. Thus, for the general iT-channel case, we provide an expression to the side distortion by approximating the 
theta series. 

Finally we will make an experimental comparison between the ^4^-i-based IA and the index assignment of 
|[T2l . which represents the state of the art in the literature. Our main focus is on the performance gain due to the 
side lattice translation. 
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A. Performance Analysis for the Two-description Case 



We first consider the description rate for encoding the Gaussian source X. Assuming high-rate quantization, the 
per-channel rate R can be shown to take the form (see 0]], 0) 

i?«ilog 2 (2^e)-log 2 (2C). (36) 

It is seen that the rate R is a function of the step size Q; if £ is known the rate is independent of the index 
assignment. The minimum rate required to transmit central indices is R c m ^ log 2 (27re) — log 2 (-^). Thus, the rate 
overhead is 2R - R c « R c - 21og 2 (2M). 

Next we study the central and side distortions, respectively. The central distortion -D(2,2) is determined by the 
central quantizer, given as 

c 2 

Note that the descriptions for the two-channel case are balanced, resulting in identical side distortions. The side 
distortion -D(2,i) can be approximated as 



D (2,l) = Y] / fx(x)(x - a (\ c )) 2 dx 
\ c eA c Jv ^ 

^ D{2 > 2) + Wl ^ (( A c-ao(A c )) 5 
A c eK.(0) 

+(A c -a!(A c )) 2 ) 



X c eV r (0) i=0 



(38) 



The part of SD costs jj J2\ ev (0) nas a simple expression 



-1 V A 2 = (l-— ) C - 
M ^ c 1 M 2j 12 

Kev T (o) 



(39) 



The part of SSD costs in (f38]> are related to J(X) by 



M-l 



£ J>(A C ) 2 = £ J(X,). (40) 

A c GK.(0) i=0 i=0 



Combining £[6]>, (|29]> and (00]) yields 

M-l M-l 

2 V " ' 2 



- «^-f). (4.) 
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Finally, inserting d37T ), d39l) and (14~T1 ) into (|38l ) produces 

M 2 C 2 

^(2,1) ~ -g^- (42) 

We now study the performance loss of the A\ -based IA. Computing the product of the side and central distortions 
from d3"7T ) and (1421) gives Z^ 2,1 ) L)( 2 > 2 ) fa This suggests that the product is unrelated to the parameter M which 
serves as a trade-off factor between side and central distortions. As £ is determined by the rate R, the product can 
be rewritten as a function of R, as given as 

The approximation (|43T > is valid for any ratio of central and side distortions as long as high-rate assumption holds. 
The work in [23] derived an approximation to the Gaussian MD rate-distortion function, which is given by 

%2)%i) ~ \2- 2R , R -»• oo. (44) 
It is seen that the gap between (l43l and (1441 is constant, which is characterized by . 

6. Performance Analysis for the K -description Case 

Consider high-rate quantization for the /^-description case. The rate R remains the same as (l36l ). Similarly, the 
central distortion D(k,k) takes the same form as Dr 2 ,2) > ie > 

D <**) « 12^- (45) 
Under the high-rate assumption, the side distortion D(k jK \ is approximated as 



1 V - 



M 

A c GK.(0) 



e E-.^- w 

v ' A c ev;(o) »=o 

The last term in (|46T > is closely related to the theta series 1 - ^ usrn § tne same analysis as in 0, the side 
distortion can be approximated as 



D (K,k) ^ ^_ 

12 



1 H -K K - 1 — -M*- 1 

A; G(5i) 



(47) 



where G(Si) is the normalized second moment of a sphere in the i dimensional space. Note the derivation from 
(|46l ) to ( 1471 involves an approximation of the theta series, which is not necessary for the two-channel case. 
Similarly to the two-channel case, some theoretical results for if-channel MD have been obtained in the past. 
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When only the individual side distortions (corresponds to k = 1) and central distortion are of primary concern in 
an MD system, the symmetric Gaussian MD rate-distortion function has been derived in [25]. The work of |fl6l 
analyzed the rate-distortion function and derived a simple approximation. Specifically, the MD rate-distortion can 
be approximated as lfT6l 

D(K,K) (D (K ,i)) K ^ *{K- lf-i K - K 2-™ R , (48) 

when R — > oo. In the next subsection for experimental comparison, we will take the approximation (l48l as a 
reference, which we denote as R-D opt. . 



C. Experimental Evaluation 

It is known that the construction of translated scalar lattices would bring the staggered gain to an MD system. 
Thus, in principle, the Ak-i -based IA would outperform the index assignment of |[T2l in the same quantization 
space. In this subsection, we evaluate the performance gain of the Ak-i -based IA over that of Ifl2l . 
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Fig. 4. A trade-off between D (K ' K ^ and D^- 1 ' for the index assignment based on the Ak-i lattice and the index assignment of [12] for 
different values of K and R = 4 (bits). The parameter M in the three experiments takes values of M = {1, 2, . . . , 7}, M = {1, 2, . . . , 10} 
and M = {1, 2, . . . , 14}, respectively. 



In order to focus on the index assignments, we consider high-rate quantization. Since we are interested in high- 
redundancy case (small M), the high-rate assumption is reasonable when the per-description rate R is large. When 
K = 2, we tested the central distortion Dr 2 ,2) m < l37b versus the side distortion D r 2j i) in (l42l) . For multi-channel case, 
we tested the central distortion D^ K K ^ in (|45T ) versus the individual side distortion Dr K ^\ in (l46l ) for K = 3, 4, 
respectively. We used the same operating rates for the index assignment of lfl2l . The experimental results are 
presented in Fig. [U one subplot for each number of descriptions. For comparison, we also plot the individual side 
distortion approximation (l42l ) or d47T ) in respective subplots. 

It is seen from the figure that when M is small, the performance gain of the A^_i-based IA over that of [12] due 
to lattice-translation is considerable. For the special case that M = 1 in the three experiments, the gain is above 2 
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dB. Further, the performance gain increases along with K. This observation suggests that for general If -description 
quantizer, it is worth to use the structure of translated lattices. 



small M. The real performance of the A^-i-based IA fluctuates closely around the derived approximation (HR- 
Approx. in Fig. @]). This suggests that in practice, one can use the approximation (l47l ) to configure the MD system 
(finding the optimal M) for a particular channel condition. 



We conclude that the proposed index assignment provides a performance gain that is significant from the point 
of view of practical MDC applications that operate with a non-vanishing redundancy among the descriptions. The 
labeling function a based on the Ak-i lattice exploits the staggered gain. It also leads to low operational complexity 
of scalar MDC as virtually the index assignment can be computed at hand. The use of the Ak-\ lattice facilitates 
analysis of the rate-distortion performance of a If -description scalar quantizer and analytic derivation of distortion 
for any description loss scenario. As a result, the obtained MDC scheme can be analytically optimized with respect 
to channel conditions enabling an instantaneous re-optimization of the scheme. 

Our results demonstrate that for the considered optimality criterion, there exist many index assignments schemes 
that are equivalent in terms of their performance. In principle any scheme from the group of the optimal index 
assignment schemes may be selected. However, certain schemes may have properties that make them particularly 
attractive. In particular, the proposed index assignment for the two-description case generates a periodic pattern of 
side quantization cells. This may lead to a low operational complexity of the resulting two-description quantizer or 
facilitate further extensions such as dithering (e.g. ifTTl ). 



Finally, we conclude from the figure that the side distortion approximation D^ K ^ in (|47T ) is accurate even for 



V. Conclusion 



Appendix A 



Proof of Lemma I3T21 



Proof: To prove the proposition, we first prove that d20])- (l2Tb take the form 



K-i r 




(49) 



3=0 



subject to 




j + z 
K 



J) 



+(2j-if+l)| =0. 



(50) 



We apply the induction argument to prove this. 
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It is obvious that d49l)-([50l) hold for z = 0. Next, we study the case that z = 1. Eq. ( f20l ) can thus be rewritten as 



K-l 



-y 

K ^ 

i=0 

1 K ~ 1 

-Y 

K ^ 



1 2 



i=0 
K-2 



K ^ 



j=0 

1 

K 

1 K ~ l 

-Y 



K(-Xi + (2(i-l)-K + l)± 



K( ■ x j+ i + (2j -K + l) 



K 



3=0 



KC(x mod(1+K _ hK) - L^-^ — -J) + {K - 1)| 



which is consistent with (|49l ). Eq. (1211 ) for z = 1 is expressed as 



K 

-y. 



i=0 
K-l 



A' 



AC-x; + (2i- A + l)| 



0. 



i=0 

Following the same derivation as for that of (l49l ). the above equation can be rewritten in the form of (l50l ). Next we 
assume the two equations hold for z = k. Similarly, we can derive the expressions for z = k + 1 based on those 
for z = k. As for z G Z~, the argument is the same. This implies that (I49l)-(r50l) hold for any zgZ. 

Note that d49t— (|50l) have the same structures as (fT4l-([T5l). By proper variable replacement, the equivalent form 
of (|49l ) is J\gc(X), where X is as given by d23l ). This completes the proof. ■ 



Appendix B 
Proof of Proposition 13.41 



Proof: The dual lattice A* K _ t is defined as EH 



K-l 



where 



i=0 



% i — K i — K 



A" ' K' K ' ' K 



with K — i components equal to i/K and i components equal to (i — K)/K, are called glue vectors. 
Let e = (1, 1, • • • , 1) T , which has K components. We rewrite the translating vector s as 



(51) 



(52) 



(53) 
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Next we define a shifting operation on a vector as 

SH([a ,--- ,a K -2,a K -i] T ) = [air ■ ■ ,a K - 1 ,a ] T . (54) 

Let sm = SH l (s) denote the vector that is obtained by performing i times the shifting operations on s. With the 
aid of d53l ). we can easily conclude that 

s [i\ =s + [i]. 

This shows that the theta series of Ak-i is identical to those of the translated lattices 1 < i < K—l. 

Thus, we have that 

©A„_>) = ^0^+4(4 (55) 
As s is a deep hole of A* K _ V this proves the proposition. ■ 
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